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Abstract 


This is the fourth semi-annual report for NAS5-97039, covering the time period July through December 
1998. The overall objective of this project is to improve the understanding of coupling processes be- 
tween atmospheric chemistry and climate. Model predictions of the future distributions of trace gases in 
the atmosphere constitute an important component of the input necessary for quantitative assessments of 
global change. We will concentrate on the changes in ozone and stratospheric sulfate aerosol, with em- 
phasis on how ozone in the lower stratosphere would respond to natural or anthropogenic changes. 

The key modeling tools for this work are the AER two-dimensional chemistry-transport model, the AER 
two-dimensional stratospheric sulfate model, and the AER three-wave interactive model with full chem- 
istry. For this six month period, we report on a modeling study of new rate constant which modify the 
NO x /NO y ratio in the lower stratosphere; sensitivity to changes in stratospheric water vapor in the future 
atmosphere; a study of N2O and CH4 observations which has allowed us to adjust diffusion in the 2-D 
CTM in order to obtain appropriate polar vortex isolation; a study of SF 6 and age of air with compari- 
sons of models and measurements; and a report on the Models and Measurements II effort. 
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(I) Evaluation of New Rates for NOy Partitioning Reactions 

Recent airborne measurements of NO and NOy obtained during the POLARIS campaign, mostly 
at high northern latitudes in summer, indicate that the NO x /NOy partitioning calculated by models for 
the lower stratosphere is systematically underpredicted. Recent measurements of reaction rate constants 
which control NO x /NOy partitioning offer a possible resolution to the model/measurement disagree- 
ment. In conjunction with David Fahey’s group at NOAA in Boulder who made the NO and NOy mea- 
surements during the POLARIS campaign, we have performed model sensitivity studies to evaluate the 
effects of these rate revisions. A paper by Gao et al. (1999) has been submitted to Geophysical Research 
Letters detailing the measurements and model calculations of NO x /NOy ratio. 

In the summer months at high latitudes, heterogeneous reactions do not play a large role in con- 
trolling the NO x /NOy ratio, Under these conditions, the ratio is largely controlled by the following reac- 
tions: 


OH + NO HN0 3 + M (Rl) 

OH + HN0 3 -> NO3 + H 2 0 (R2) 

HNO3 + hv — > OH + N0 2 (R3) 

The rates for reactions Rl and R2 reported in JPL-97 (DeMore et al., 1997) have recently been ques- 
tioned. Donahue et al. (1997) first suggested that the rate for reaction Rl reported in JPL-97 may be too 
high. Dransfield et al. (1998) and Brown et al. (1998) reported new values for the rate of reaction Rl for 
stratospheric conditions which are about 20% lower than JPL-97 values. Brown et al. (1998) also re- 
ported values of the rate constant for R2 which are higher than the JPL-97 recommendation. Using a 
photochemical steady state box model, Gao et al. (1999) determined that agreement between models and 
measurements of NO x /NOy is much improved with the Brown et al. (1998) rates for Rl and R2 relative 
to the JPL-97 rates (model/measurement ratio increases from 0.62 to 0.84). 

The Brown et al. (1998) rates for Rl and R2 have been incorporated into the AER 2-D chemical 
transport model. Figure la shows the change in calculated ozone for June of 2015 with subsonic aircraft 
included due to the modified rates relative to the JPL-97 rates. Ozone reductions of up to 5% are seen 
between 20 and 40 km altitude, while ozone increases are seen in the troposphere. Changes in ozone 
column (Figure lb) exceed 1% only in northern polar regions. The percentage of ozone loss due to NO x 
increases by 2-8% in the lower stratosphere, while that due to HO x and C10 x decreases. This change in 
the relative importance of the ozone loss cycles modifies the sensitivity of ozone to perturbations. 

A new rate constant has also been measured for the reaction of N0 2 with O: 

N0 2 +0->N0+0 2 . (R4) 

This reaction does not affect the NO x /NOy partitioning significantly, but is a loss mechanism for odd 
oxygen. The new rate (A. R. Ravishankara, personal communication, 1998) is about 30% higher than 
the JPL-97 value at 200 K. Figure 2 shows the calculated changes in ozone due to inclusion of this 
modified reaction, along with Rl and R2, in the AER 2-D model. Local ozone is reduced by up to 13% 
at 30 km in the high latitudes, up to 6% at the equator at 30 km, and increased in the troposphere. Col- 
umn ozone is reduced by 1-7% relative to a model with the JPL-97 rates. 

We have calculated the ozone change due to HSCT emissions in 2015 (with 3.0 ppbv of Cly) and 
2050 (with 2.0 ppbv of Cly) under two different emission scenarios. Both scenarios include 500 HSCTs 
operating at Mach 2.4 at 18-20 km altitude with an emission index (El) for NO x of 5, but scenario 1 in- 
cludes no sulfur emissions, while scenario 2 includes sulfur emissions assuming EI(SO 2 )=0.4 and 50% 
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conversion of the emitted sulfur to small particles in the aircraft wake. Table 1.1 summarizes the calcu- 
lated annual average ozone column changes at 60°N. Without sulfur emissions, the Brown et al. rates 
for R1 and R2 lead to a slight increase in ozone depletion due to HSCT because of the increased impor- 
tance of the NO x cycle for ozone loss. With sulfur emissions, the new rates have the opposite effect. 
Ozone depletion due to emitted sulfur (which increases the stratospheric aerosol burden) is a function of 
the HO x and C10 x cycles, which are reduced under the new chemistry. Including new rates for reac- 
tions Rl, R2, and R4 leads to even more ozone depletion due to HSCT for scenario 1 and even less 
ozone depletion for scenario 2. With Rl, R2, and R4 updated, column ozone depletion at 60°N in 2015 
is the same with or without sulfur emissions, and in 2050, column ozone depletion is greater without sul- 
fur emissions than with sulfur emissions. We will be examining the the effect of modifications in Rl, 
R2, and R4 on the sensitivity of ozone response to chlorine within the next few months. 


Table 1: Annual average column ozone change at 60°N due to a fleet of 500 HSCTs operating at Mach 
2.4 with EI(NO x )=5 in the year 2015 or 2050. 


HSCT 

EI(SO?) 

with JPL-97 rates 

Rl, R2 updated 

Rl, R2, R4 updated 

Scenario 


2015 

2050 

2015 

2050 

2015 

2050 

1 

0 

-0.5% 

-0.5% 

-0.6% 

-0.6% 

-0.8 

-0.7 

2 

0.4 

-1.3% 

-0.9% 

-1.0% 

-0.5% 

-0.8 

-0.4 


(II) Evaluation of Modified Stratospheric Water Vapor 

For the 1998 NASA report on the impact of High Speed Civil Transport (HSCT) aircraft on the 
atmosphere, AER was asked to perform calculations with modified concentrations of stratospheric water 
vapor. Our standard treatment of H 2 O fixes concentrations in the troposphere based on relative humid- 
ity and calculates H 2 O in the stratosphere, using an imposed boundary condition at the tropopause. 
Stratospheric water vapor normally has mixing ratios of 3-6 ppmv. For these sensitivity studies, the 
tropopause boundary condition for H 2 O was decreased by 2 ppmv, increased by 2 ppmv, and increased 
by 4 ppmv, resulting in commiserate changes in H 2 O throughout the entire stratosphere. Increases in 
H 2 O result in increases in HO x , which increases ozone loss, particularly in the lower and upper strato- 
sphere where HO x forms the dominant ozone loss mechanism. Changes to background ozone concen- 
trations under these three perturbations are shown in Figure 3 for a 2015 atmosphere with subsonic air- 
craft. Decreases of 2 ppmv in stratospheric H 2 0 result in increases in ozone column of 1-4%, with 
larger changes in polar springtime. Increases of 2 ppmv in stratospheric H 2 O results in decreases in 
ozone column of 1-3%, with larger changes in polar springtime. Increases of 4 ppmv in stratospheric 
H 2 0 result in ozone column decreases of 2-7% outside the polar springtime regions, up to 12% in the 
northern polar spring and 15% in the southern polar spring. Local ozone changes in the upper strato- 
sphere due to modified water are +12%, -8%, and -14% for H 2 O decreases of 2 ppmv, increases of 2 
ppmv, and increases of 4 ppmv, respectively. Changes in ozone in the mid-latitude lower stratosphere 
are as much as +8%, -8%, and -18% for H 2 O decreases of 2 ppmv, increases of 2 ppmv, and increases 
of 4 ppmv, respectively. 

The ozone response to HSCT is sensitivity to the background stratospheric water vapor concen- 
tration. At mid to high latitudes, especially in summer and fall, HO x represents the dominant loss 
mechanism for ozone in the lower stratosphere. Decreases in H 2 O (and HO x ) leads to an increase in the 
importance of the NO x loss cycle in this region, resulting in a greater ozone response to HSCT. How- 
ever, in the polar winter and spingtime, increases in H 2 Q concentration make the triggering of PSC 
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events by HSCT-emitted N0 X more likely. Figure 4 shows calculated changes in column ozone due to 
HSCT for the scenario with 500 HSCTs operating at Mach 2.4 with EI(NO x )=5 and no sulfur emissions 
in 2015. For the case with H 2 0 reduced by 2 ppmv, maximum ozone depletion in the northern hemi- 
sphere occurs in late summer and fall; PSC events are not triggered. With the normal H 2 0 concentra- 
tion, maximum depletion of 1 .0% occurs in spring, but depletions of 0.6% are seen in summer also. 

With H 2 0 increased by 2 ppmv, summer ozone depletion is only -0.2%. With H 2 0 increased by 4 
ppmv, springtime ozone depletion of up to 2% is calculated, and summertime shows incomplete recov- 
ery from the spring depletion. The sensitivity of ozone to increases in chlorine under these conditions of 
modified stratospheric water will be examined in the coming months. 


(Ill) POLARIS Campaign and Analysis of 2-D CTM Transport 

In conjunction with our participation in the POLARIS campaign, we have been analyzing the 
transport rates of the AER 2-D model and comparing modeled and derived quantities with available ob- 
servations. Transport affects the ozone distribution directly by regulating the concentration of ozone and 
indirectly by influencing ozone production and loss through the concentrations of source and reservoir 
gases. Recent analyses using calculated SFg concentrations along with the derived mean age of air have 
been performed to test four different versions of transport within the AER 2-D CTM. Three of these 
have different horizontal mixing rates to simulate horizontal transport ranging from a weak tropical bar- 
rier to a strong tropical barrier. Recent analyses from the Models and Measurements II intercomparison 
exercise have indicated that many of the atmospheric models have circulations that are too vigorous 
(Hall et al. 1998). In response to this finding, one of the transport schemes scales the stream function by 
0.6, creating slower velocities in both the horizontal and vertical directions. A comparative analysis of 
these models and in situ measurements was performed and general impacts on the HSCT assessment 
calculations were examined. 

Using Geller et al.’s (1997) mean global tropospheric growth rate (equation 1) of SF$ for the 
years 1987 through 1996, based on a wide range of tropospheric measurements, 

[SF 6 ](t) = 3.4361 + 0.2376(t - 1996) + 0.0049(t - 1996) 2 (1) 

and Volk et al.’s (1997) derived equation which directly relates the SF6 (pptv) concentration to a mean 
age of air (equation 2), 


[SF 6 ](x,t) = [SF 6 ](t-r-8D + 2cA(r-80. (2) 

SF 6 concentrations and the mean age of air for the four different transport versions were calculated. 

Here, [SF 6 ](x,t) is the SF^ concentration in pptv at point (x,t) in space and time, 8T is the tropospheric 
lag time (0.8 years), T is defined as the mean age of air in years (first moment of the Green’s function), 
c is 0.0049 pptv/yr 2 and A (1.25 years) parameterizes an appropriate value of the width of the age spec- 
trum by using a weighted average of two GCM outputs. Measurements of SFg and N 2 0 concentration 
were used from the Stratospheric Tracers of Atmospheric Transport (STRAT) (Elkins, et al., 1996), the 
Airborne Southern Hemisphere Ozone Experiment/Measurements for Assessing the Effects of Strato- 
spheric Aircraft (ASHOE/MAESA) (Tuck, et al. 1997), and the Photochemistry of Ozone Loss in the 
Arctic Region In Summer (POLARIS) campaigns to compare the different transport versions of the 
model. 
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Figure 5 shows the calculated 1996 annual average age of air in years as a function of the annual 
average N 2 0 concentration (ppbv) for each of the versions of the AER 2-D CTM with the in situ mea- 
surements overlaid. Three regions of the models are represented: the southern mid-latitudes, the tropics 
and the northern mid-latitudes. In all versions of the model, the northern mid-latitude region calculates 
the oldest ages of air. The youngest ages of air are seen in the tropical region. This is directly at- 
tributed to the large scale transport, with a weaker hemispheric motion moving the SF^ concentrations 
into the northern hemisphere and large scale flow from the troposphere entering the stratosphere via the 
tropical tropopause. 

Figure 6 illustrates the annual average 1996 calculated mean age of air as a function of latitude 
and N 2 0 concentration with the in situ measurements overlaid. Younger air is seen at high tropical alti- 
tudes in the LEAKY-PIPE (panel b) and the PIPE (panel c) versions of the model, indicative of a stron- 
ger tropical upwelling in these versions of the model. The LEAKY-PIPE and PIPE show similar struc- 
ture in the age contours, with the PIPE version calculating the youngest ages overall. The NO-PIPE 
(panel a ) and the SLOW (panel d) versions of the model resemble one another qualitatively but the 
SLOW model calculates ages of air that are much older than the NO-PIPE version. The SLOW transport 
version (panel d) closely resembles the measurements, both qualitatively and quantitatively, calculating 
the flattest contours in the lower altitudes and the oldest ages of air in the higher latitudes and altitudes 

Figure 7 shows the change in calculated column ozone as a function of latitude and season for a 
500 aircraft fleet with emission index of 5, flying in a 2015 atmosphere, with background aerosols. In- 
spection of Figure 7 shows the largest ozone depletion due to the HSCT perturbation, in the SLOW 
(panel d) version of the model. Total column ozone change in the SLOW version of the model for the 
northern mid-latitudes ranges from -0.5% to -0.8% loss in contrast to the value of — 0.3% in the NO- 
PIPE (panel a) version of the model. The larger ozone loss in the NH is directly due to the change in the 
transport quantities, where a stronger barrier coupled with a slower stream function prevents much of the 
aircraft exhaust to circumvent the globe. The LEAKY-PIPE, and PIPE ozone responses, shown in pan- 
els b and c respectively, fall between the SLOW and NO-PIPE responses. 

Although the SLOW version of the model gives reasonable agreement with the in situ measure- 
ments of SF 6 concentration, its calculated ages, and various other tracer species, it is not necessarily a 
consistently better transport scheme in all regions of the atmosphere. Table 2 shows calculated global 
lifetimes in years for selected tracers from all versions of the AER 2-D model. Also shown in Table 2 
are calculated ages, based on in situ measurements, from Volk et al. (1997) Comparison of the global 
ages derived from in situ measurements with the four different transport versions of our model give 
clues to transport issues. All of the versions of the model give different ages, with the PIPE giving the 
youngest and the SLOW giving the oldest ages. Notably, the SLOW version calculates ages that are er- 
roneously long for all species not dependent on OH chemistry. 

Based on Figure 5, which shows age of air vs N 2 0 for the four versions of the AER 2-D CTM, 
the question arises why the SLOW model does not match Volk’s calculated lifetime for N 2 0. Careful 
analysis indicates that the age of air, in the AER model, does not share the same mixing surfaces with 
other long-lived species. This does not indicate to what extent the real atmosphere is in steady state and 
whether true SFg concentrations share the mixing surfaces of other long-lived species in the atmosphere. 
It does however, merit careful consideration when trying to extrapolate theory to the real atmosphere. 
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Table 2: Calculated global lifetimes (years) for selected species from the AER 2-D models. Lifetimes 
calculated by Volk et al. (1997), based on correlations with age of air, are included for comparison. 


SPECIES 

SLOW 

NO-PIPE 

LEAKY-PIPE 

PIPE 

VOLK et al. 

n 2 o 

171.0 

131.2 

112.3 

102.1 

124+/- 49 

C2C13F3 

122.1 

94.2 

78.7 

71.8 

89 +/- 35 

CBrClF 2 

52.2 

38.1 

33.8 

32.3 

20 +/- 9 

CFCI3 

74.0 

56.4 

48.4 

45.3 

41 +/- 12 

CF 2 C1 2 

147.1 

112.4 

95.0 

85.8 

77 +/- 26 

CC1 4 

63.0 

48.0 

41.6 

39.3 

32 +/- 1 1 

CH3CCI3 

7.0 

6.1 

6.0 

5.8 

- 

ch 4 

12.7 

11.3 

11.2 

11.0 

- 


(IV) Model Adjustment of Mixing Rates across the Polar Barriers 

The results of the Models and Measurements II experiments showed that there are significant de- 
fects in the transport parameters currently used in most 2-D chemistry transport models. One of these 
defects is the lack of a clear polar vortex: Large amounts of observations and measurements show that 
the fast isentropic mixing created by the planetary wave breaking cannot penetrate into the polar vortex 
and the concentrations of long-lived species, e.g. N2O and CH4, have strong gradients rapidly across the 
polar vortex well. This feature is not seen in the results of the AER 2-D CTM (Figure 8a). The isen- 
tropic mixing is represented by the eddy diffusion coefficient, Kyy in the 2-D CTM. Smaller values of 
Kyy should be assigned in the polar barrier to simulate the polar vortex. The measurements of long- 
lived species by instruments aboard of UARS with nearly global coverage and high precision provide a 
useful database for adjusting Kyy across the polar vortex. Recently, Randel et al. (1998) processed the 
UARS/MLS and UARS/CLAES CH4 data to create the altitude-equivalent latitude distributions of CH4 
by averaging the data along the potential vorticity (PV) contours to avoid smoothing the gradients of 
CH4 distribution across the vortex edge. These new 2-D distributions of CH4 reveal the polar barriers, 
which greatly reduce the mixing between the mid-latitudes and polar region (Figure 8b). 

During the past 6 months, we have used these CH 4 distributions to adjust the mixing across the 
polar vortex in the AER 2-D CTM. We define a cost function to measure the difference between the 
gradients of the CH 4 on the polar vortex edge from the model simulations and from Randel’s data. By 
adjusting the Kyy to minimizing the cost function calculated using a 2-D CTM model with parameter- 
ized CH 4 chemistry, we derive values for Kyy near the polar vortex. The values of Kyy, which produce 
the minimum for the cost function have significantly improved the model simulations of CH 4 across the 
polar barrier (Figure 8c). These adjusted values of Kyy for both north and south polar vortices were 
used in the 2-D CTM with detailed photochemistry. The new version of the full 2-D model improves 
the simulations of other long-lived species in the polar region. Unlike in the subtropical barrier case, the 
adjustment of Kyy in the polar barrier only has localized effects on the age of air calculation and the 
simulated atmospneric effects of the High Speed Civil Transport (HSCT). The exchange time scale be- 
tween the polar vortex and the middle latitudes is calculated using the same method as for the subtropi- 
cal barrier (Shia et al, 1998). This study was presented at the 1998 AGU Fall Meeting and is reproduced 
here as Appendix B. 
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(V) Activity Associated with the Models and Measurements Workshop II 

After the March 1998 meeting, drafts of different sections were collected from various authors. 
Malcolm Ko, Charles Jackman, Alan Plumb and Karen Sage met at NASA GSFC on October 6 for a one 
day meeting to finalize the structure of the report. The draft report was sent to NASA Langley in De- 
cember for copy editing. 
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Figure 1 : Calculated percentage change in (a) local ozone in June, and (b) total column ozone due to 
modifying the rates for reactions R1 and R2 from JPL-97 values to those reported by Brown et al. 
(1998). Results shown are for 2015 with subsonic aircraft emission included. 
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Figure 2: Calculated percentage change in (a) local ozone in June, and (b) total column ozone due to 
modifying the rates for reactions R 1 , R2, and R4 from JPL-97 values to those reported by Brown et al. 
(1998) and Ravishankara (personal communication, 1998). Results shown are for 2015 with subsonic 
aircraft emission included. 
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Figure 3: Calculated changes in June ozone (left column) and total ozone (right column) due to changes 
in stratospheric H 2 O relative to the standard treatment in the AER 2-D model. Panels (a) and (b) have 
stratospheric H 2 0 reduced by 2 ppmv, panels (c) and (d) have stratospheric H 2 0 increased by 2 ppmv, 
and panels (e) and (f) have stratospheric H 2 0 increased by 4 ppmv. All calculations represent the 2015 
atmosphere with a subsonic aircraft fleet. 
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Figure 4: Calculated changes in total column ozone due to an HSCT fleet of 500 aircraft operating at 
Mach 2.4 with EI(NO x )=5 in 2015 with nonvolcanic aerosol surface area SAO. Panel (a) is with strato- 
spheric H 2 O decreased by 2 ppmv, panel (b) with stratospheric H 2 O at normal levels, panel (c) with 
stratospheric H 2 O increased by 2 ppmv, and panel (d) with stratospheric H 2 O increased by 4 ppmv. 
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Figure 5: Calculated 1996 annual average age of air in years for all latitudes and altitudes as a function 
of N 2 O concentration (ppbv). Insitu measurements are represented by symbols. Profiles from 35°N, the 
Equator and 35°S are represented with various line styles. The NO-PIPE version is displayed in panel 
(a), the LEAKY -PIPE version in panel (b), the PIPE in panel (c) and the SLOW in panel (d). 
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Figure 6: Annually averaged 1996 calculated age of air (years) from SF 6 concentrations for the four 
AER 2-D CTMs as a function of latitude and N 2 0. Filled circles represent the corresponding in situ 
measurements. Model contours are overlaid. The NO-PIPE version is displayed in panel (a), the 
LEAKY-PIPE version in panel (b), the PIPE in panel (c) and the SLOW in panel (d). The contour in- 
terval is 0.25 years. 
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Figure 7: Calculated change in column ozone (%) due to an HSCT perturbation scenario of EI(NO x )=5 
with 500 aircraft for a 2015 atmosphere with background aerosols. The NO-PIPE version is displayed in 
panel (a), the LEAKY-PIPE version in panel (b), the PIPE in panel (c) and the SLOW in panel (d). The 
contour interval is 0.1%. 
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Figure 8: Latitude-altitude distribution of CH4 mixing ratio (ppmv) in the northern hemisphere for Feb- 
ruary (a) as simulated by the AER 2-D CTM, (b) from UARS observations (Randel et al., 1998), and (c) 
as simulated by the AER 2-D CTM with polar barrier. 
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Abstract. Extensive airborne measurements of the reactive ni- 
trogen reservoir (NO y ) and its component nitric oxide (NO) 
have been made in the lower stratosphere. Box model simula- 
tions that are constrained by observations of radical and long- 
lived species and which include heterogeneous chemistry sys- 
tematically underpredict the NO x (= NO + N0 2 ) to NO y ratio. 
The model agreement is substantially improved if newly meas- 
ured rate coefficients for the OH + HN0 3 and OH + N0 2 reac- 
tions are used. When included in 2-D models, the new rate co- 
efficients significantly increase the calculated ozone loss due 
to NO x and modestly change the calculated ozone abundances 
in the lower stratosphere. Ozone changes associated with the 
emissions of a fleet of supersonic aircraft are also altered. 


Introduction 

Understanding the mechanisms controlling the abundance 
of NO y and its partitioning into component species is an essen- 
tial requirement for understanding the stratospheric ozone (0 3 ) 
layer. Reactions involving NO x form a catalytic O3 destruc- 
tion cycle and also moderate 0 3 loss due to other cycles in- 
volving reactive hydrogen (HO x ) and halogens (CiO x -BrO x ) 
[cf. Wennberg et aL, 1994]. The partitioning of the NO y reser- 
voir between NO x and other component species involves gas- 
phase and heterogeneous reactions as well as photolytic proc- 
esses (see Figure 1) [cf. Gao et aL, 1997]. We present here ex- 
tensive new measurements of NO and NO y obtained in the 
high latitude lower stratosphere during summer. This data set 
provides a unique opportunity to test our understanding of the 
gas-phase chemistry linking NO x and nitric acid (HNO3), 
which is generally the most abundant NO y species. Because the 
continuous daylight present at summer high latitudes limits 
the heterogeneous production of HNO3 by N 2 05 hydrolysis, 
gas-phase reactions primarily control the balance between NO x 
and NO y . Outside summer polar regions, the N 2 05 hydrolysis 
reaction occurring on stratospheric sulfate aerosols is a more 
important sink of NO x [cf. Fahey et aL, 1993], particularly in 
the lower stratosphere during winter when heterogeneous 
pathways account for most of the HNO3 production [cf. Gao et 
aL , 1997]. 

Observations and box model description 

The observational data set used here was acquired with in- 
struments on board the NASA ER-2 high altitude aircraft dur- 
ing the Photochemistry of Ozone Loss in the Arctic Region In 
Summer (POLARIS) mission. The mission used deployment 
sites at latitudes of 37°N, 65°N, and 24°N. In the present 
study, data from 23 flights are used to examine NO x /NO y at 
high latitudes (> 60°N) in late spring, summer, and early fall 
periods, and in the tropics in early fall. 
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Measurements of NO and NO y are made with a three- 
channel chemiluminescence detector | Gao et al., 1997]. Be- 
cause N0 2 measurements are not available for all flights, N0 2 
values inferred from the steady state relationship (N0 2 *) [ Gao 
et al ., 1997] are used throughout for consistency. The agree- 
ment between average N0 2 * and observed N0 2 values during 
POLARIS is within 10% (L. A . Del Negro et al., Comparison of 
modeled and observed values of N0 2 and JNO2 during the 
POLARIS mission, submitted to J. Geophys. Res., 1998). The 
uncertainty of the NO x /NO y measurements is estimated to be 
±20%. Measured aerosol surface area (SA) densities varied be- 
tween 0.5 and 1.5 pm 2 cnr 3 in the data set used here. 

A photochemical steady state box model that includes only 
the processes shown in Figure 1 is used to predict NO x /NO y in 
a sampled air parcel. The reaction set used in the model is a 
subset of the more comprehensive set used by Salawitch et al. 
[1994]. The model is constrained by observed values of NO, 
OH, CIO, 03, SA, pressure and temperature, and column O3 
above the aircraft [see references in Gao et al., 1997 for instru- 
ment details]. Because NO, OH, and CIO concentrations ap- 
proach zero at high solar zenith angles (SZAs), data gathered at 
SZA > 85° are not used in this work. The diurnal dependence 
of the OH radical was determined empirically and normalized 
here to values observed along the flight track [Wenrtberg et al., 
1994; T. F. Hanisco et al ., unpublished data, 1999]. Modeled 
CI0N0 2 values agree with POLARIS in situ observations to 
within 20% (R. M. Stimpfle et al The coupling of C10N0 2 , 
CIO, and N0 2 in the lower stratosphere from in-situ observa- 
tions using the NASA ER-2 aircraft, submitted to J. Geophys. 
Res., 1998). Model BrO values are estimated by calculating 
the steady state partitioning of the Br y reservoir as estimated 
from organic bromine observations [Wamsley et al., 1998]. 
Unless otherwise noted, gas-phase rate coefficients, absorption 
cross sections, and reactive uptake coefficients on aerosol are 
taken from NASA JPL-97 recommendations [DeMore et al., 
1997]. Photolysis rates are calculated using a radiation scatter- 
ing model [Salawitch et al., 1994] which includes effects of 
overhead O3, albedo and cloud heights. Reactive uptake coef- 
ficients of 0.1 and 0.8 are used for N 2 C>5 and BrON0 2 , respec- 
tively. The model does not include C10N0 2 hydrolysis since 
the reactive lifetime in sampled air parcels exceeds 100 days 
and thus this process has a negligible effect on NO x /NO y 
[Robinson et al., 1997]. With input parameters averaged or in- 
terpolated to 100 s intervals, the model is run to a diurnal 
steady state to yield a value of NO y consistent with measured 
NO and rate parameters used in the model. The NO x /NO y de- 
duced in this way is then compared with that inferred from 
N0 2 * and the measured values of NO and NO y . 

Results 

The comparison of observed and modeled NO x /NO y shown 
in Figure 2A, includes all stratospheric data from altitudes be- 
tween 17 and 20 km (50 mb < pressure < 80 mb) with tempera- 
tures greater than 204K. The approximately 200,000 s of data 
are distributed in latitude between 0° and 90°N (< 40°: 5%; 
40° - 60°: 17%; > 60°: 78%). In high latitude summer, NO x 
reaches 3 parts per billion by volume (ppbv) with NO x /NO y in 
the range of 0.2 to 0.3. Although similar NO x /NO y ratios can 
also be found in the tropical lower stratosphere, the associated 
NO x values are typically much less than 3 ppbv. 

Most of the steady state model values of NO x /NO y fall sig- 
nificantly below observed values. The linear fit (forced 
through the origin) to the model/data regression has a slope of 
0.65. The regression remains virtually unchanged if JNO2 is al- 
tered to account for the small difference between N0 2 * and ob- 
served N0 2 values because JNO2 changes affect N0 2 values on 
both axes in Figure 2A. 



The diurnal steady state assumption is unlikely to be the 
cause of the observation-model discrepancy in Figure 2A. 
Back trajectories of the sampled air parcels were calculated us- 
ing National Center for Environmental Prediction (NCEP) 
analyzed winds and temperatures. For each parcel, the accumu- 
lated latitude change over the previous 10 days was calculated 
using a weighting function inversely related to the time before 
aircraft sampling. When those parcels with an effective latitude 
change of 5° or more (~ 70% of the data) are excluded from the 
data set, the regression fit to the remaining data is not signifi- 
cantly changed. 

The model parameterization of N 2 0 5 hydrolysis is also un- 
likely to be the cause of the NO x /NO y discrepancy. Using 
NCEP back trajectories, sampled air parcels that had recently 
experienced nearly continuous solar illumination were identi- 
fied. The acceptance criterion was an SZA of less than 93° for 
90% of the preceding 5 day period. As solar illumination be- 
comes continuous, the diurnal conversion of NO x to HNO3 
through N2O5 formation and hydrolysis becomes negligible 
[Farman et al, f 1985; Briihl et al., 1998] and hence the con- 
tribution of N2O5 hydrolysis to NO x /NO y also becomes negli- 
gible. The separately averaged data sets in Figure 2A, show 
that, overall, these selected data points are nearly indistin- 
guishable as a group from those that experience interrupted so- 
lar illumination. In contrast to the diminished role of N2O5 
hydrolysis, the contribution of BrON0 2 hydrolysis to reduc- 
ing NO x /NO y is greatest for the continuously illuminated par- 
cels [Randeniya et al., 1997]. However, BrON0 2 hydrolysis 
reduces NO x /NO y by a maximum of only 5 - 10%. Thus, the 
hydrolysis reactions of N2O5, C10N0 2 , and BrON0 2 are not 
effective enough individually or collectively to cause the sys- 
tematic discrepancy shown in Figure 2A. 

In air parcels in which the heterogeneous hydrolysis reac- 
tions in Figure 1 do not play a strong role, NO x /NO y is primar- 
ily controlled by the gas-phase reactions: 


OH + N0 2 + M — HN0 3 + M 

0) 

OH + HNO3 -* NO3 + H 2 0 

(2) 

HNO3 + hv — OH + N0 2 

(3) 


where M is 0 2 or N 2 . Only a few measurements of the rate coef- 
ficient for (1) and (2) (kj and k 2 ) are available at lower strato- 
spheric temperatures (220 ± 20K) and pressures (50 - 150 mb) 
[DeMore et al., 1997]. Donahue et al. [1997] first suggested 
that JPL-97 values for kj may be too high for stratospheric 
conditions. Recently, Dransfield et al. [1999] and Brown et 
al. [1999] reported laboratory measurements of for strato- 
spheric conditions that are about 20 - 30% lower than JPL-97 
values but consistent with the few previous measurements 
available below 298K. Brown et al. [1999] also show that k 2 
in the lower stratosphere is higher (up to 50%) than JPL-97 
values. 

Incorporating the new values of kj and k 2 of Brown et al. 
[1999] in the model significantly improves the overall agree- 
ment with the observations (Figure 2B). The linear fit (forced 
through the origin) with the new model results has a slope of 
0.90. In a similar analysis using balloon observations of 
NO x /NO y at ER-2 altitudes and above, agreement to within 
10% is achieved for a 35% reduction of kj from JPL-97 values 
[Osterman et al., 1999]. The Brown et al. values of k| are not 
expected to substantially affect the Gao et al. [1997] results for 
winter polar NO y partitioning since N2O5 hydrolysis is gener- 
ally faster than reaction (1) under the conditions sampled. In 
contrast, an increase in k 2 will lead to a calculated increase in 
NO x /NO y during winter. 

An increase of the HNO3 photolysis rate in (3) would also 
improve the observation-model comparisons in Figure 2. In 
the absence of changes to kj and k 2 , an increase of approxi- 


3 



matcly 70% is required to yield a regression slope of 0,9. Such 
a large increase exceeds the reported uncertainties in the under- 
lying measurements of the temperature dependent HN0 3 ab- 
sorption cross section for lower stratospheric conditions 
[ Burkholder et al., 1993] or in the uncertainties associated 
with the radiation field calculations [ Gao et al ., 1997; 
Salawitch et al., 1994|. 

Implications 

The 2-D dynamical-chemical-radiative model of Garcia 
and Solomon (GS) [1994] and the AER chemistry-transport 
model [ Weisenstein et al., 1996] were used to examine the ef- 
fects of the Brown et al. rates on O3 and the principal ozone 
catalytic loss cycles at high latitudes in summer. The results 
shown in Figure 3 include altitudes between 15 and 40 km in 
the stratospheric ozone layer in July at 59°N to illustrate the 
effect of the new k 1 and k 2 values over a range of model tem- 
peratures and other conditions. The important features in Fig- 
ure 3 are (i) that large changes in the rates that control 
NO x /NO y yield much smaller changes in the ozone distribu- 
tion in the lower stratosphere due to the corresponding changes 
in HO x , C10 x , and BrO x abundances, and (ii) that the new rates 
result in a ‘cross-over’ near 20 km from an increased to a de- 
creased value of the total O3 destruction rate. Throughout 
most of the 15 - 40 km region, the magnitude of the NO x - 
catalyzed destruction increases (up to 10% near 25 km when 
expressed as a percentage of the new total O3 destruction rate) 
as a direct result of higher NO x /NO y values while the contribu- 
tions from the BrO x -C10 x and HO x catalytic cycles decrease 
because of moderation by NO x . Above 25 km, NO x /NO y ap- 
proaches unity and the sensitivity of the 0 3 loss cycles to the 
rate coefficient changes approaches zero. At 25 km, the net in- 
crease in 0 3 destruction rates reaches a maximum. Because 
changes in k 2 are largest at low temperatures and affect the 
abundances of both NO x and HO x , the sensitivity of the BrO x - 
C!O x and HO x cycles to increases in NO x /NO x is largest below 
25 km and results in a net decrease in the total 0 3 loss rate. 
Changes in the vertical distribution of O3 due to the new rates 
are small (< 5%) throughout the 15 - 40 km region, consistent 
with the calculated changes in the O3 loss rate. Brown et al. 
[1999] show that lower-latitude O3 changes are also small (< 
2%) in the GS model using the new kj and k 2 values and an in- 
creased value for the O + NC^ rate coefficient. 

The increased contribution of the NO x cycle to 0 3 loss rates 
at mid and high latitudes modifies how emissions from a pro- 
posed High Speed Civil Transport (HSCT) stratospheric aircraft 
fleet are expected to change 0 3 in the lower stratosphere 
[Stolarski et al., 1995]. Whether the Brown et al. rates will 
lead to a larger or smaller depletion for a particular HSCT sce- 
nario depends sensitively on the relative emissions of NO x and 
sulfur from an HSCT engine. The sulfur emission index (El) 
and exhaust plume processes determine the SA changes of the 
background sulfate aerosol. As an example. Table 1 shows how 
the Brown et al. kj and k 2 values affect O3 column amounts at 
59°N in summer for two key HSCT fleet scenarios in a 2015 
and 2050 atmosphere. With no sulfur emissions, the changes i n 
the rate coefficients of kj and k 2 result in slightly larger O3 
depletion because the role of NO x in catalytic O3 loss is greater 
than its role in moderating the HO x and halogen cycles. When 
sulfur is also emitted, the formation of sulfate particles in the 
plume and their accumulation in the atmosphere causes an in- 
crease in background SA densities. Larger SA densities reduce 
NO x values through N 2 C>5 hydrolysis (for both the NO x emit- 
ted by the aircraft and NO x in the background atmosphere), re- 
sulting in a decrease in the 0 3 removal by the NO x loss cycle 
compared to the atmosphere without aircraft. The NO x reduc- 
tion due to increased SA densities in the atmosphere increases 



the HO x and active halogen concentrations, resulting in in- 
creases for the corresponding O 3 removal cycles. The effects of 
increased NO x concentrations due to aircraft emission must be 
considered in conjunction with the effects of increased SA den- 
sities. The use of the Brown et al. rates in this case results in less 
O 3 depletion because more effective moderation of the HO x 
and halogen loss cycles by aircraft -emitted NO x partially miti- 
gates the effect of aerosol changes. 


Concluding remarks 

The in situ aircraft observations presented here represent a 
large new data set with which to examine stratospheric 
NO x /NOy over a wide range of latitude and season in the lower 
stratosphere. Representative box model calculations systemati- 
cally underpredict the observed ratio by about 40%. Because 
the discrepancy remains in air parcels with continuos solar ex- 
posure, the parameterization of the N 2 O 5 hydrolysis reaction is 
unlikely to be the cause. The use of recent laboratory results 
for the OH + N0 2 and OH + HN0 3 reaction rate coefficients 
from Brown et al. [1999] significantly improves the compari- 
son of measured and modeled NO x /NO y for most sampled air 
parcels. The effect of the rate changes and associated NO x in- 
creases on O 3 concentrations calculated in 2-D models is small 
(< 5 %) throughout most of the lower stratosphere because of 
the interdependence of the NO x , Br0 x -ClO x , and HO x cata- 
lytic loss cycles. 
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Table 1. Column O 3 change (%) at 59°N in July 3 


EI(S0 2 ) (g/kg fuel) 

JPL 97 rates 

Brown et al. rates for (1) and (2) 

0 

-0.5/-0.4 

-0.6/ -0.6 

0.4 b 

-1 .3/0.9 

-1.0/ -0.5 


a AER model results for the scenario of 500 HSCT aircraft operating 
at Mach 2.4 near 18 km cruise altitude with an emission index (El) of 
(NO x ) of 5 g/kg fuel. Values are noted for both 2015 and 2050 models 
years as ‘2015/2050’. Similar values are calculated for average hemi- 
spheric 0 3 column changes. 

b With 50% conversion to small sulfate particles in the plume. 
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Figure 1, Schematic of the reaction pathways between the 
principal N0 y component species in the lower stratosphere. 
Photolysis reactions are indicated by hv. ‘Sulfate aerosol’ de- 
notes heterogeneous reactions on background aerosol particles. 

Figure 2. Comparison of modeled and measured values of 
NO x /NO y from POLARIS. The symbols are averages of data as 
grouped by increasing NO x /NO y value. Open symbols are av- 
erages of 80 points, each of which represented 100 s of observa- 
tional data. The solid symbols are separate averages (84 points 
of 100 s data) for those POLARIS observations for which the 
sampled air parcel experienced near-continuous solar illumina- 
tion for the preceding 5 days. The horizontal and vertical bars 
on each symbol represent the l-o sample standard deviation 
within each group of measured and modeled values, respec- 
tively. Results are shown using the model with JPL recom- 
mended rate coefficient values (panel A) and with the Brown et 
al. [1999] rate coefficient expressions for (1) and (2) (panel B). 
The solid and dashed lines represent the 1:1 and 1:2 regression 
slopes, respectively. 

Figure 3. GS model results vs altitude for percentage changes 
in ozone, its total destruction rate, and the contributions of the 
NO x , HO x , and halogen catalytic cycles for average July 1990 
conditions at 59°N with background aerosol conditions. The 
changes are those that result from using the new Brown et al. 
[1999] values for kj and k 2 . The percentage changes of the 
catalytic cycles are shown with respect to the new total ozone 
loss rate. The percentage changes between the new and old val- 
ues of the individual loss cycles are larger. For example, at 20 
km the NO x , HO x , and halogen cycle changes are +23%, -10%, 
and -22%, respectively. Nearly identical results are obtained 
with the AER model. The representativeness of the GS model 
was checked by comparing the average observations at 65 °N 
with the nearest model grid point (64°N). Noon model values 
of NO x /NO y are somewhat lower (< 25% ) in the 15 - 20 km re- 
gion at high latitudes. The longer-lived NO y ,Cl y ,and O3 spe- 
cies are within 10 - 20% of the average observed values be- 
tween 15 and 20 km. 


Figure 1. Schematic of the reaction pathways between the principal NO y component species in the lower stratosphere. Photolysis 
reactions are indicated by hv. ‘Sulfate aerosol’ denotes heterogeneous reactions on background aerosol particles. 

Figure 2. Comparison of modeled and measured values of NO x /NO y from POLARIS. The symbols are averages of data as grouped 
by increasing NO x /NO y value. Open symbols are averages of 80 points, each of which represented 100 sof observational data. 
The solid symbols are separate averages (84 points of 100 s data) for those POLARIS observations for which the sampled air parcel 
experienced near-continuous solar illumination for the preceding 5 days. The horizontal and vertical bars on each symbol repre- 
sent the l-o sample standard deviation within each group of measured and modeled values, respectively. Results are shown using 
the model with JPL recommended rate coefficient values (panel A) and with the Brown et al. [1999] rate coefficient expressions 
for (1) and (2) (panel B). The solid and dashed lines represent the 1:1 and 1:2 regression slopes, respectively. 

Figure 3. GS model results vs altitude for percentage changes in ozone, its total destruction rate, and the contributions of the 
NO x . HO x , and halogen catalytic cycles for average July 1990 conditions at 59°N with background aerosol conditions. The 
changes are those that result from using the new Brown et al. [1999] values for k( and k 2 . The percentage changes of the catalytic 
cycles are shown with respect to the new total ozone loss rate. The percentage changes between the new and old values of the indi- 
vidual loss cycles are larger. For example, at 20 km the NO x , HO x , and halogen cycle changes are +23%, -10%, and -22%, respec- 
tively. Nearly identical results are obtained with the AER model. The representativeness of the GS model was checked by com- 
paring the average observations at 65°N with the nearest model grid point (64°N). Noon model values of NO x /NO y are somewhat 
lower (< 25% ) in the 15 - 20 km region at high latitudes. The longer-lived NO y , Cl y , and O3 species are within 10 - 20% of the 
average observed values between 15 and 20 km. 
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